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^S) ! Abstract. We calculate the optical Hall conductivity for mono- and bilayer graphene 

within the Kubo formalism and demonstrate that the ratio of the jump in the 
Hall conductivity in bilayer graphene is twice the corresponding value of monolayer 
f/2 I graphene. Furthermore, we find a characteristic logarithmic singularity as the 

frequency approaches the gap energy. The optical Hall conductivity is almost 
unaffected by thermal fluctuations and disorder in monolayer graphene. Only for 
bilayer graphene disorder has a stronger effect on transport properties. 
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'^ I Monolayer graphene (MLG) is a monoatomic sheet of carbon atoms arranged in a 

Q I honeycomb lattice with unique transport properties. This is a consequence of the two- 

dimensional nature of the material and due to the band structure which consists of 
two separate bands touching one another at isolated nodal points. In the vicinity of 
>■ I the nodal points quasiparticles exhibit a linear spectrum. The main difference between 

MLG and bilayer graphene (BLG) is that the low-energy excitations of the latter have 
^ I a quadratic spectrum [H |2]. For the longitudinal conductivity this difference causes a 

factor of 2 for the DC conductivity [3] and also for the optical conductivity |1]. This 
O ' leads to the question how a change of the low-energy spectrum affects the quantum HaU 

' properties of few layer graphene. 

An intriguing phenomenon in graphene is the quantum Hall effect (QHE) which 
Sp ' ^^^ already observed in the first experiments on graphene [5]. It exhibits a rather 

;h ' unexpected anomalous behavior. In contrast to the QHE in a two-dimensional electron 

gas of a semiconductor, the Hall plateaux appear antisymmetrically around zero carrier 
density p. Additionally, the magnitude of the Hall conductivity of the first plateau is for 
BLG twice the corresponding value of MLG [1]. This doubling of the Hall conductivity 
is similar to the above mentioned effect found for the longitudinal conductivity. An open 
question is whether or not the doubling appears also for the optical Hall conductivity. 
This shall be studied in the rest of this paper using the Kubo formalism. 

It is widely accepted that the QHE occurs in semimetals as a consequence of a gap 
opening and broken time reversal symmetry. Usually, this is achieved by applying a 
magnetic field perpendicular to the 2D plane, whereas the gap opening alone can be 
obtained by hydrogenation of MLG [6] or by a double gate in the case of BLG [2]. 
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Since the broken time reversal symmetry implies a broken symmetry between the 
two nodal points in MLG or BLG, we can study the QHE by assuming the absence of 
intervalley scattering and focusing only on one nodal point. Such a situation can be 
created by applying a periodic magnetic field [7] or a spin texture [8]. 

In the following we calculate the Hall conductivity for low-energy quasiparticles 
near a nodal point with a uniform gap and a power-law spectrum with an integer 
exponent. Then the low-energy Hamiltonian, describing electrons in graphene-like 
systems with a uniform gap A = 2771, reads in Fourier representation as 

where n = 1 should be associated with MLG and n = 2 with BLG. For general n the 
eigenvalues of the Hamiltonian are 



Ei = {-iyE, E = -fVm^ + P" , (2) 

where Z = (/ = 1) refer to the upper (lower) band, respectively. The corresponding 
eigenvectors are 




V^.^(r) = ^^^( ^^r™ lexp(zk-r). (3) 

The band parameter 7 for MLG is 7 = f i^' [9] and for BLG 7 = Vp/'ji [10], where vp is 
the Fermi velocity and 71 ~ OAeV is the interlayer coupling constant. To simplify the 
notation we drop indices for spin- and valley degeneracy and put 7 equal to unity. 

The Hall conductivity can be calculated as the off-diagonal element of the Kubo 
conductivity tensor [1]: 

, , {Ei\j,\E,){E,\j,\Ei) 



' flV 



lim - / y^ 



oh ^^ El -El, 



J{Ei,-EF)-f{Ei-EF) d^k 

Ei-Ev+u-ia (27r)2' ^ '' 

where Ep represents the Fermi energy, f{E) = 1/(1 + exp(/3£')) the Fermi-Dirac 
distribution at the inverse temperature /3 and u the frequency of the external field. 
The current operator 

jf, = ie[H,rf,] (5) 

has vanishing diagonal elements. This is due to the fact that for cu 7^ the excitations 
El, El, satisfy \Ei — Ei,\ = u, and it is also reflected by vanishing intraband matrix 
elements. First we calculate current matrix elements defined in (jll). Due to rotational 
symmetry of the model, the use of polar coordinates is more convenient. Since the 
angular variable enters only the current matrix elements, the corresponding integration 
can be carried out separately. Then we obtain for the intraband matrix elements 

>'2tt 

{±E\j^\±E){±E\jy\±E)dip = 

4e^ r^2fc4n-. cos(,.)sin(^) ^^ ^^ 

. n 4:E'^ 



Optical Hall conductivity of gapped graphene 3 

The nonvanishing interband contribution of the matrix elements reads 

J {±E\j,,\TE){TE\jy\±E)dcp = ±2e't7Tn'e--'-. (7) 

In contrast to the calculation of the longitudinal optical conductivity in |1], the current 
matrix elements are now imaginary. In order to obtain the real part a'^^ of 0"^,^ we have 
to evaluate a Cauchy principal value integral 



a 



^ f ST^ o 2 • 2 7 2n-2 "^ 

— / > Ze mn k — 

f{Ev - Ep) - j{Ei - Ep) kdk 
El -Ev+uj (27r)^ 
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Substituting the k integration by an integration over the energy E = \^rn? + W^, the 
corresponding Jacobian becomes 

J=f-V' = -^ (9) 

All powers of k in the integrand cancel each other, such that the real part of the Hall 
conductivity reduces to the simple expression 

, ne'm f^ f(-E-Ep)-f{E-Ep) 

The imaginary part is given by 

ne'^m [^ f{-E - Ep) - f{E - Ep) 



a" 



xy 



TXl 



h ./i^i 2E 



x[5{2E + u)-5{2E~oj)]dE. (11 

In the limits T — t- and Ep — > we obtain 

2ira + oj 



, n e^ m 



2m — uj 



(12) 



„ ne^ m 



<y'iy = -—-e{uj-2m), (13) 

which reduces in the DC limit w — ;■ to 

2 
Tl 6 

(^'xy = sgn(m) - — , a'^y = 0. (14) 

This result was also found by Semenoff for n = 1 [11]. Hence, the DC Hall conductivity 
is a nonzero constant proportional to the spectral power exponent n. We would like to 
emphasize once again that the time reversal symmetry breaking is crucial for observing 
the QHE. If we had considered a second nodal point (i. e. with time reversal symmetry), 
the Hall conductivities from both cones would cancel each other. 

For large u the optical Hall conductivity decays like u~^. Remarkable is the 
singularity at a; = 2m in ( IT2l) . It appears when the external frequency u is equal 
to the gap of the electronic system. The latter means that there is an additional current 
contribution to the transport due to the excitation of states in the upper band. This 
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is a kind of resonance that disappears for higher frequencies. Moreover, instead of 
the singularity the longitudinal optical conductivity only jumps to a finite value at 
u = 2m |1]. The singularity of the optical Hall conductivity could be used to determine 
the gap experimentally. 

In figure [1] we show a plot of expression ( TTOl) as a function of u and different 
temperatures. The Hall conductivity scales with n. This is in agreement with QHE 
experiments, where the plateau of BLG around zero carrier density is twice that of 
MLG [11 [5]. The temperature dependence is controlled by the energy scale of the 
corresponding system. In MLG, the relevant energy scale is the hopping parameter t 
~ 2.8ey, which corresponds to a temperature ~ 32.5 x IQ^K. The relevant scale of BLG 
is even higher, since it is oc t^/71. Therefore, expression fITOl) is insensitive over a wide 
range of temperatures T ^ Tp, as shown in figure [H In particular, the gap singularity 
of the optical Hall conductivity survives. For very high temperatures, however, the Hall 
conductivity is reduced and goes eventually to zero. This is shown for fixed frequency 
uj in figure El 
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Figure 1. Temperature dependence of the Hall conductivity. 



Disorder: We have seen that thermal fluctuations have almost no effect on the 
singularity at w = 2m. Since graphene is also subject to disorder effects (ripples, 
impurities, etc.), we study their influence on the singularity in the following. Returning 
to the Kubo formula (|1]), we can rewrite the conductivity as (cf. [13] ) 

{Tr[j^6{Hd,s-e')j,6{Hd,s-e)]) 
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(15) 



e- e 

where (...) represents the disorder average. The latter can be approximated in the self- 
consistent Born approximation by replacing the Hamiltonian Hdis by {Hdis) + if] [H]. 
The average Hamiltonian {Hdis) is the same as the Hamiltonian in ([T]) and 77 is the 
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Figure 2. Optical Hall conductivity near the singularity as a function of the inverse 
temperature /3. 



scattering rate caused by the disorder. This imphes that we have to replace in ( ITSjl the 
Dirac delta functions by 



S{H, 



dis 



S^{H 



271 



[{H - e + iT])-^ - {H - e - iT])-^] . (16) 



In MLG the scattering rate is ?7 oc exp{7r/g) [15], where g is the variance of the random 
gap. In BLG, on the other hand, rj is proportional to the variance g, which implies that 
the influence of disorder on BLG is much stronger. Realistic fluctuations in graphene 
are less than a tenth of the hopping rate g ^ 0.1. This implies that a realistic scattering 
rate would be r/ ^ 2 x 10^^^ (MLG) and r/ ^ 0.1 (BLG). 

Now the trace can be expressed again in the diagonal representation as 

Tr[jMH-e')jJ,{H-e)] = 

krlk 

{E,\j,\E,){E,\j^\Ei) 6r,{Ei - e')6,{E, - e) ^^ • (17) 

The angular integration can be carried out (cf. ([7])) and transforming the /c-integral to 
an energy integral gives 

W-Ep)-fp{e-Ep) 1 



e^nm 
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xy 



27ih J|„| J_^ J_^ e-e' e-e' + UJ 

X [5r,{E - e')6r,i-E - e) - 6r,{-E - e')5r,{E - e)] dede'dE. (18) 

We evaluated expression (TT8l) for T = and /i = numerically for several values of rj. 
The results are depicted for the real and imaginary part of the optical Hall conductivity 
in figures [3] and HI One can see that the Hall plateaux remains nearly constant for all r] 
under consideration, whereas the singularity is broadened. 

Conclusions: In this work we have studied the optical Hall conductivity of gapped 
graphene within a low-energy approximation. Our calculations indicate that the DC 
Hall conductivity in MLG and BLG only depends on the sign of the mass term and on 
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Figure 3. Real Part of the optical Hall conductivity for different values of the 
scattering rate r]. 
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Figure 4. Imaginary part of the optical Hall conductivity for different values of the 
scattering rate ?/. 



the exponent of the low-energy spectrum and it reproduces the experimentally observed 
factor of 2. Interesting is that the optical Hall conductivity is quite insensitive to 
thermal fluctuations over a wide range of temperatures. The effect of the curvature of 
the low-energy spectrum is also surprisingly simple: The optical Hall conductivity of 
MLG and BLG differ by a factor 2, as it was also found for the longitudinal optical 
conductivity. The same factor 2 appears in the visual transparency of MLG and 
BLG [12]. Interestingly, there is a logarithmic singularity in the optical Hall conductivity 
when the frequency u of the external AC field becomes equal to the gap of the electronic 
system. The appearance of the singularity in our calculations is related to particle- 
hole excitations and can be interpreted as a resonance phenomenon. Although thermal 
fluctuations have no effect on this singularity, disorder may soften it in the case of 
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BLG, where the scattering rate can be large. In case of MLG, where the scattering 
rate is very small, the singularity of the optical Hall conductivity is almost unaffected. 
Consequently, the singularity could be used to determine the gap in MLG by measuring 
the optical Hall conductivity. 
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